A sample of Seyfert-2 galaxies with ultra-luminous galaxy-wide 

NLRs - Quasar light echos? * 
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ABSTRACT 

We report the discovery of Seyfert-2 galaxies in SDSS-DR8 with galaxy-wide, ultra-luminous 
narrow-line regions (NLRs) at redshifts z = 0.2 — 0.6. With a space density of 4.4 Gpc -3 
at z ~ 0.3, these "Green Beans" (GBs) are amongst the rarest objects in the Universe. We are 
witnessing an exceptional and/or short-lived phenomenon in the life cycle of AGN. The main focus 
of this paper is on a detailed analysis of the GB prototype galaxy J2240-0927 (z = 0.326). Its 
NLR extends over 26 x 44 kpc and is surrounded by an extended narrow-line region (ENLR). With 
a total [Olll]A5008 luminosity of (5.7=1=0.9) x 10 43 erg s -1 , this is one of the most luminous NLR 
known around any type-2 galaxy. Using VLT/XSHOOTER we show that the NLR is powered 
by an AGN, and we derive resolved extinction, density and ionization maps. Gas kinematics is 
disturbed on a global scale, and high velocity outflows are absent or faint. This NLR is unlike any 
other NLR or extended emission line region (EELR) known. Spectroscopy with Gemini/GMOS 
reveals extended, high luminosity [O ill] emission also in other GBs. WISE 24/im luminosities are 
5 — 50 times lower than predicted by the [O ill] fluxes, suggesting that the NLRs reflect earlier, 
very active quasar states that have strongly subsided in less than a galaxies' light crossing time. 
These light echos, or ionization echos, are about 100 times more luminous than any other such 
echo known to date. X-ray data are needed for photo- ionization modeling and to verify the light 
echos. 

Subject headings: Galaxies: active — Galaxies: evolution — Galaxies: Seyfert 
1. Introduction 
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Ministerio da Ciencia, Tecnologia e Inovacao (Brazil) and 



NLRs are a common sight in active galaxies 
and explained by the unified AGN model (see e.g. 
Antonucci|[l993| |Bianchi et al.||2Q12[ ). UV/X-ray 
emission from the central black hole engine photo- 
ionizes the ISM over large distances, but may be 
shielded by neutral gas or dust absorbers. The 
prototype NLR consists of two ionization cones 
(e.g. in NGC 4151 and NGC 5252), whic h may 
fragment into several individual clouds (|Evans| 



et al. 1993 Tsvetanov et al. 1996). In Seyfert 



2 galaxies NLR sizes are a few hundred pc, and 



Ministerio de Ciencia, Tecnologia e Innovacion Productiva 
(Argentina). 
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can extend up to 5 kpc (Bennert et al. 2006a). 



Emission line FWHMs are on the order of several 
hundred km s _1 and are often richly structured. 

Extended narrow-line regions (ENLRs) were 
first described by Unger et al.| ( |1987| ) and have 
been found around many AGN since (Durret 



1989). They extend over 20 kpc or more and have 
much lower luminosities than NLRs. Dynamical 
FWHMs are low (< 50 km s _1 ), whereas high 
excitation levels reveal AGN as ionizing sources. 
ENLRs usually follow the galactic rotation. 
Extended emission line regions (EELRs 



Fos- 



bury 1989) are found mostly around radio- loud 



QSOs and have similar sizes as ENLRs, but very 
different characteristics. Fu fc Stockton| ( |2009 ) 
find compact clouds with low line widths (< 100 
km s _1 ) but moving with high velocities (~ 500 
km s _1 ). Dynamically chaotic structures are com- 
mon, and morphological links to the host galaxies 
absent. EELRs are likely the result of a merger, 
kick-starting the QSO engine that ionizes the gas 
and blasts it into the outer surroundings. 



Lintott et al. (2009) find a cloud of ionized gas 



20 kpc outside IC 2497. Lacking an apparent AGN 
in the host galaxy, it is interpreted as a light echo 
from a very active earlier AGN phase. Indeed ra- 



dio observations by Rampadarath et al. (2010) re- 
veal hidden AGN features. |Keel et al.| fl2012a| ) 
present 154 galaxies with similar detached clouds 
at redshifts z < 0.1, finding evidence for short 
periods of high AGN activity. During these peri- 
ods the luminosity can change by up to 4 orders of 



magnitude, as has been shown by|Schawinski et al. 



(2010) and Keel et al. (2012b) for IC 2497. How- 



ever, time-scales inferred by Keel et al. (2012a 
are too short to be explained by current accretion 
models. Such light echos therefore provide new 
insights into the onset and shutdown processes of 
QSO activity on scales of a galaxy's light crossing 
time (rsj 10 4 — 10 5 years). 

Galaxy formation models predict AGN-driven 
large outflows needed to explain properties of 
the interstellar medium and massive galaxies (e.g. 
Springel et al.|[2005 ). NLRs on galaxy scales have 
been found in powerful radio galaxies (HzRGs, 
Nesvadba et al.|2008[), and ultra- lum inous infrared 
galaxies (ULIRGs, |Harrison et al. 2012| ) at high 
redshifts (z = 2 - 3). Line widths of 700 - 1400 
km s _1 and large offsets with respect to the sys- 
temic redshift emphasize the outflow character in 



the ULIRGs with highest [O ill] luminosity. 

|Cardamone et aL (2009) present a new class 
of emission- line galaxies, dubbed "Green Peas" 
(GPs), referring to their compact size and color 
in gri images. In a sample of 112 GPs they find 
80 star- forming galaxies, 9 Seyfert-1, 10 Seyfert- 
2 and 13 transition objects (showing both AGN 
and star formation features). GPs with AGN-like 
emission characteristics have so far not been stud- 
ied beyond the initial identification and analysis in 



Cardamone et al. (2009), and differences to NLRs 



in other Seyfert galaxies have not been found. 

It is now established that supermassive black 
holes (SMBHs) reside in the centers of massive 
galaxies, which must have hosted an AGN in their 
past. Consequently, mergers of galaxies must also 
lead to the coalescence of SMBHs. Indeed nu- 
merous wide-separation binary AGN are known, 



however systems with kpc separations (Shen et al. 



2011 Comerford et al. 2012|) or less (Rodriguez 



et al. 2006 Fabbiano et al. 2011) are rare and 



have to be verified carefully. Simulations of the 
pre-coalescence state of SMBH mergers (such as 
Hopkins Quataert|[2010l |Khan et"aT||2gl2l fVanl 



Wassenhove et al.||2012 ) do not predict extraordi- 



nary NLR properties. 

In CFHT/Megaprime data we serendipitously 
discovered J224024. 1-092748 (hereafter: J2240), a 
peculiar galaxy at z = 0.326 with GP colors. Due 
to its large angular extent of 4" x 1" J2240 was 
not included in the GP sample of | Cardamone et aT 



(2009). In this paper we show that it is fundamen- 



tally different from GPs. We have identified two 



dozen similar objects in SDSS-DR8 (Aihara et al. 
2011) with redshifts z = 0.1 — 0.7, after remov- 
ing 95% spurious detections. We refer to these 
galaxies with large and ultra-luminous NLRs as 
"Green Beans" (GBs). GBs display a previously 
unknown phenomenon in the life of AGN, and we 
explore different formation mechanisms. 

The paper is structured as follows. In Sect. 2 
we present observational data for J2240, followed 
by our 2D spectral analysis method in Sect. 3. 
We discuss the physical properties obtained from 
this analysis in Sect. 4, after which we present our 
sample of similar galaxies extracted from SDSS. 
We discuss our findings and summarize in Sect. 6. 
Throughout this paper we assume a flat stan- 
dard cosmology with Q m = 0.27, = 0.73 and 
Hq = 70/^70 km s _1 Mpc _1 . The relation be- 
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Fig. 1. — Left: CFHT/Megaprime gri color image of J2240. The photometric redshift of the fainter reddish 
galaxy to the left of the nucleus categorizes it as a background source, whereas the galaxy below with the 
prominent spiral arm has the same spectroscopic redshift as J2240. The dashed line shows the orientation of 
the VLT/XSHOOTER slit, running (from the upper left to the lower right) through the Disk and the nucleus, 
and then through the Cloud. Right: Complex appearance of the [Om]A5008 line in the VLT/XSHOOTER 
spectrum (log-scaled, continuum subtracted). Some of the features discussed in this paper are marked. The 
dashed line indicates the spatial position of the nucleus. 



tween physical and angular scales at z = 0.326 
is l'/O = 4.76 kpc. Error bars represent the la 
confidence level. All wavelengths stated are vac- 
uum wavelengths. 

2. Observations 
2.1. CFHT imaging 

Deep ugriz CFHT/Megaprime data (Table [TJ 
were obtained by us through Opticon proposal 
2008BO01 at 2008-09-20 in excellent conditions 
to study a supercluster of galaxies at z = 0.45. 
Further details about these data can be found in 
Schirmer et all (120111). J2240 at z = 0.326 is a 



serendipitous discovery in these images, extend- 
ing over 7" x 4". With a pixel scale of 0"186 and 
Of! 7 image seeing the galaxy is well resolved. It 
has irregular morphology and peculiar colors sim- 
ilar to that of a GP (see Fig.[T]), with the exception 
of very low ix-band flux (u — r — 4.06 mag). K- 
corrections and luminosities were calculated using 
kcorrect v4.2 flBlanton fc Roweis||2007|). Most 



noteworthy are a bright non-stellar nucleus in a 
disk-like body reminiscent of a spiral galaxy, and 
a half-detached cloud about 12 kpc south-west of 
the nucleus. The cloud extends over 8 x 18 kpc. 



We hereafter refer to these two components sim- 
ply as the 'Disk' and the 'Cloud', without implying 
any properties about their physical shape. Spec- 
troscopic observations reveal a much more com- 
plex picture, as can be seen in the right panel of 
Fig. [I] We refer to the central ±4 kpc encompass- 
ing the nucleus as the 'Center'. 

Figure [I] shows a fainter reddish galaxy with a 
photometric redshift of z = 0.37 ± 0.06, projected 
2'/0 east of the nucleus of J2240. A second galaxy 
with a prominent spiral arm is located 4V 5 to the 
south-east, and has the same spectroscopic red- 
shift as J2240 (z = 0.326). 

2.2. Spectroscopy 

2.2.1. VLT/FORS2 

The redshift of J2240 and its classification 
as a Seyfert-2 galaxy were first secured using a 
VLT/FORS2 ( |Appenzeller et al.||1998"| spectrum 
(ESO DDT proposal 286.A-5027). Our data were 
reduced using a custom pipeline. In short, all spec- 
tra were overscan corrected, debiased, flat-fielded, 
rectified and sky subtracted. Flux-calibration and 
correction for telluric absorption were based on 
the central white dwarf of the Helix Nebula (NGC 
7293) and the atmospheric transmission model for 
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Table 1: Characteristics of the CFHT/Megaprime 
data and J2240 photometry. The limiting magni- 
tudes represent the 5<r completeness limit for non- 
stellar sources. Luminosities are dereddened. 
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Paranal ( |Patat et al.||2011 ). A technical summary 
of the characteristics of our spectroscopic data sets 
is given in Table [2] 

2. 2. 2. VLT/XSHO O TER 

Using VLT/XSHOOTER ( |Vernet et al.1|2QlT| ), 
we obtained a spectrum (proposal 287.B-5008) of 
the full 300 — 2500nm range, with a slit posi- 
tion angle of 43.5 degrees (same as for FORS2) 
aligned along the major axis of the extended emis- 
sion in Fig. [I] Spectral (spatial) resolutions are 
5-20 (1.5) times higher than with FORS2. The 
data were pre-processed using the XSHOOTER 
pipeline ( Modigliani et al. 2010). Flux calibra- 



tion (using the spectrophotometric standard Feige 
110), continuum subtraction and further data pro- 
cessing were done using custom-made software. 
Since the slit width for the flux standard was much 
wider than the one for the science target, and the 
standard was calibrated with a different flat field, 
the absolute calibration is only known to within a 
constant factor. We determine the latter from the 
FORS2 spectrum, over which we have full con- 
trol and consistent calibrators. After correction, 
FORS2 and XSHOOTER fluxes are indistinguish- 
able within their noise in the common wavelength 
area. The full XSHOOTER spectrum is shown in 
Appendix [Aj Fig. [12) 

2.2.3. Gemini/ GMOS 

We used Gemini-South/GMOS (proposal GS- 
2012A-Q-83) to create a redshift survey of further 
GB candidates (Sect. |5|. Given that our targets 
are bright and our only aim was to detect bright 
emission lines, this campaign was designed as a 
backup program to be executed in bad seeing and 



cirrus. Most data were taken in the second half 
of the 2012A semester. At the time of writing 
observations are still ongoing. 

2.2.4- Correction for galactic extinction 

J2240 (GLAT = -54°) lies in an area of faint 
galactic cirrus. However, the particular line of 
sight is unaffected and has E(B — V) = 0.059 mag 
and R = 3.3 flSchlegel et al.|1998[ ). The extinction 
models for R = 3.3 and the 'standard' R = 3.1 are 
indistinguishable in the wavelength range relevant 



for this paper ( |Fitzpatrick 1999). Hence we use 
the R = 3.1 extinction curve flCardelli et aT|l989[ 



Osterbrock & Ferland 2006). Correction factors 



are between 1.09 for (redshifted) [S n] A6718,33 and 
1.25 for [Nev]A3427. 

3. 2D spectroscopic analysis method 

J2240 is a dynamically complex system, and 
well resolved in our data. A proper description 
requires models of the NLR's 3D structure and 
knowledge about the radiation field, serving as in- 
put for photo-ionization codes. This is beyond the 
scope of this discovery paper, as we are currently 
lacking the X-ray data to characterize the ioniz- 
ing spectrum, as well as IFU observations of the 
entire object. Nevertheless, applying 2D emission 
line diagnostics we can obtain information going 
beyond single global values for e.g. extinction or 
density. To do so, we have to bring all emission 
lines to a common reference system. 

The emission lines are complex (see e.g. the 
right panel of Fig. [I]). We have therefore decided 
to use a parameter-free approach and work with 
the line images directly, as compared to model- 
ing with multiple Gaussians. For example, a su- 
perposition of up to four Gaussians of different 
widths is needed to describe any spatial line scan 
through the Cloud, which fragments into several 
smaller condensations. Parts of the lines are in ad- 
dition asymmetric, generally making a Gaussian 
parametrization a bad choice in this case. 

3.1. Theoretical [OIII] and [Nil] doublet 
line ratios 

The theoretical value of the [O lll]A4960,5008 
intensity ratio has long been underestimated com- 
pared to observations (e.g. 2.89 by |Galavis et al.| 
19971). |Leisy fc Dennefeld| ([1996]) found 3.00±0.08, 
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Table 2: Characteristics of the spectroscopic data. GMOS data were taken throughout the 2012A semester. 
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consistent with high-S/N measurements of other 
nebulae, and encouraged improved theoretical cal- 
culations to be carried out. 

|Storey fc Zeippen ([2000) calculate precise line 
ratios for [O ill] A4960, 5008 and homologous cases 
(such as [Nn]A6550,86, [Nein]A3869, 3968 or 
[Ne v]A3347, 3427) using relativistic corrections. 
The intensity ratios of these lines are independent 
of temperatures and densities in typical astro- 
physical nebulae, as the two lines come from a 
common upper level. Only if the lines become op- 
tically thick will the intensity ratio change. This 
is unlikely, however, given their low absorption co- 
efficients. If the density was high enough to cause 
significant optical depth, collisional de-excitation 
would suppress these lines in the spectrum (Pe- 
ter J. Storey, priv. comm.). The ratio of the 
transition probabilities between [Olll]A5008 and 
[Olll]A4960 is found to be 3.013, which translates 
into a flux ratio of 2.984 by multiplying with the 
wavelength ratio 4960/5008. This value has been 
confirm ed (2.993 ± 0.014) by pimitrijevic et al. 
fl2007[) based on 34 SDSS AGN spectra. Like- 



wise, for [Nll]A6586 and [Nll]A6550 the flux ratio 
evaluates to 3.05. 

We use these theoretical values to deblend the 
Ha- [Nil] line complex, and to verify our spectral 
analysis method with the [Olll] doublet. 

3.2. Re-projection of lines 

As a first step we subtract the continuum, 
which we linearly interpolate across an emission 
line based on the adjacent ~ 100 pixels on both 
sides. The resulting continuum is smoothed in 
dispersion direction with a 100 pixel wide median 
kernel, and subtracted from the original spectrum. 
We do not correct line fluxes for potential stellar 
absorption, as the equivalent widths are high (sev- 



eral 100 A to more than 1000A). 

In the second step we extract images centered 
on each line of interest, encompassing their full 
spatial extent and 2. 3 A in width. These images 
are then stretched by a factor of Aii ne /AH/3 to cor- 
rect for the wavelength dependence of the velocity 
broadening. We choose H/3 as the reference frame, 
as it is in the middle between the blue [Oil] and 
red [S n] lines. The re-projection onto H/3 also 
includes a 2-fold binning along the dispersion di- 
rection to increase the S/N. We retain an effective 
spectral resolution of R ~ 3800. 

Thirdly, we have to register the line images such 
that they overlap precisely. Distortion correction 
by the XSHOOTER pipeline left no measurable 
spatial offsets of the continuum recorded in the 
UVB and VIS arms. Registering the lines in dis- 
persion direction is more difficult, as their individ- 
ual appearances are different. Fortunately, almost 
all lines exhibit traces of the ENLR with low line 
widths, which we use as a reference mark. Re- 
maining lines are registered using other common 
features. In this way registration in dispersion di- 
rection is accurate to ~ 1 pixel, well within the 
oversampled spectral resolution. 

Lastly, we smooth with a 2 pixel wide Gaus- 
sian kernel, cut off at a radius of 2 pixel. Prior 
to smoothing, the highest and lowest pixel in the 
aperture are rejected. In this way we suppress 
spurious noise features, which can get strong when 
calculating ratios with small denominators. 

Error maps from the XSHOOTER pipeline are 
treated analogously and fully propagated. 

3.3. Deblending close line pairs 

Velocity broadening in J2240 leads to some 
overlap between Ha and [N il] A6550,86, and within 
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-■(!■) [Nil] A6550 [NII]A6586 





(2) Subtracting H/3from (1), extracting [Nil] A6586 



(3) Subtracting [Nil] from (1) --> Ha 




(4) Subtracting (3) from (1) --> [Nil] 




Fig. 2. — Deblending the Ha- [Nil] line complex 
(panel 1). Panel 2: A rescaled and smoothed ver- 
sion of H/3 is subtracted, and [Nll]A6586 extracted 
within a closed polygon line. Panel 3: [Nil] gets 
smoothed and subtracted from the original line 
complex, leaving Ha isolated. Panel 4: Ha gets 
smoothed and subtracted from the original line 
complex, leaving the [Nil] doublet. 



the [Sll]A6718,33 doublet. While the bright cores 
of these lines are well separated, their fainter wings 
overlap. We deblend them as follows. 

For the Ha- [Nil] complex (Fig. |5J panel 1) we 
have to perform several steps. First, we smooth 
H/3 with a 2 pixel wide Gaussian and subtract it 
from Ha after multiplication with 3.0. The latter 
is near the average Case B hydrogen ratio (3.08) 
in Seyfert-2 AGN flGaskell fc Ferland|[l984l ). In 
this way we remove about 90 — 95% of the Ha 
contamination in the [Nil] lines (Fig. |2j panel 2), 
thus preparing a first guess of the [Nil] doublet. 
Some residual Ha is still present, occurring due to 
variable extinction in J2240 (Sect. [XT] ). 

Next we extract [Nll]A6586 within a closed 
polygon line, then smooth and subtract it from 
the original Ha- [Nil] complex. We also divide 
it by 3.05 (Sect. |3.1[ ), shift it to the position of 
[Nll]A6550, and subtract it once more. The [Nil] 



doublet has thus been removed and we have pre- 
pared a clean Ha line (Fig. |5J panel 3). Lastly, 
we smooth the Ha image and subtract it from the 
original line complex, leaving the final [N n] pair. 

For [S ll]A6718, 33 we use a different approach, 
as this doublet is reasonably separated with little 
overlap. We define non-overlapping closed polygo- 
nal lines of identical shape around both lines. Pix- 
els within each polygon line are smoothed with a 
small Gaussian kernel and subtracted, leaving the 
other line in the doublet isolated. 

Note that our 2D analysis technique used in the 
following is confined to regions with high S/N. Ar- 
eas where the [S il] doublet or Ha and [N n] overlap 
are already faint and mostly excluded. 

3.4. [OIII] sanity check 



As shown in Sect.pU] the [O lll]A4960,5008 line 
ratio is independent over a large range of tempera- 
tures and densities. Both lines are strong and un- 
contaminated, and thus serve as a good test case 
to validate our 2D analysis method. While a con- 
stant offset indicates a problem with flux calibra- 
tion, systematic variations along spatial or spec- 
tral directions hint at bad registration of the two 
lines, or problems in the data processing or extinc- 
tion correction. 

Figure [3] shows that in our data the [O ill] ra- 
tio is constant within errors, thus validating our 
2D analysis method. Utilizing the area where 
jA5008/jA4960 uncertainty is less than 0.03, we ini- 
tially estimated a line ratio of 3.169 ± 0.012 (in- 
cluding full extinction correction, see Sect. 4.1). 



This is inconsistent with the theoretical value, and 
amongst the highest values measured by |DimitrT 
jevic et al. (2007). Retracing our processing steps 
we identified an unfortunate coincidence causing 
this discrepancy. The [Olll]A4960 line is red- 
shifted to 6577A, just redwards of Ha. Being 
a white dwarf, the spectrophotometric standard 
Feige 110 shows a ~ 20 A wide Ha absorption 
line in the XSHOOTER spectrum. The refer- 



ence flux values tabulated by Oke (1990), how- 



ever, are based on a low resolution spectrum (~ 
5 A pixel -1 ), resulting in a measured line width of 
^ 40A. As the [Om] line is narrow (1.5A) com- 
pared to the absorption line, we can easily deter- 
mine a correction factor of 1.075 ± 0.015. The 
[Om] intensity ratio then becomes 2.950 ± 0.015. 
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Fig. 3.— The [Om] jASOOs/j^geo ratio (theoret- 
ical value: 2.98) is constant within errors. The 
black ellipse shows the spatial and spectral scale 
beyond which data points become statistically in- 
dependent (same for all other similar plots). The 
dashed line indicates the spatial location of high- 
est continuum flux (coincident with the nucleus). 
The origin of the velocity axis coincides with the 



location of highest electron density (see Sect. 4.2). 



The same effect is present in our FORS2 data (for 
which a different white dwarf served as flux cali- 
brator), yielding a corrected ratio of 3.01 ± 0.04. 
No other emission lines are affected. 

4. Results from the spectral analysis 

We use bright emission lines to determine ex- 
tinction, electron density and temperature, abun- 
dance and ionization. These diagnostics are well 
understood for various ionizing spectra. However, 
as the X-ray spectrum of J2240 is not available, 
uncertainties do exist in some of the parameters 
derived. We address this where applicable. 

4.1. Internal extinction 

4.1.1. Method 

The extinction of a line with intrinsic inten- 
sity I\o is given as I\ = I\oe~ c f( x \ where C is 
a constant, and /(A) parametrizes the extinction 
curve of the interstellar medium, and depends on 
the chemical and physical properties of the dust 
grains. Using the observed and intrinsic (before 
attenuation by dust) intensity ratio of two emis- 



E(B -V) = -2.b\og 



I x /I v \ f(B) - f(V) 



ho/ J /(A) -/(i/) 



(1) 

We use the Ha/H/3 ratio and the standard 
R = 3.1 model flCardelli et~aL]|1989| ) to calcu- 
late the reddening in J2240. This requires knowl- 
edge of the intrinsic value of Ha/H/3, which equals 
~ 2.87 for Case B conditions in a typical Hll re- 
gion. In the harder radiation field of an AGN col- 
lisional excitation of Ha; becomes important, in- 
creasing the ratio to ~ 3.08. However, it depends 
strongly on Lx/L opt and the slope of the X-ray 
spectrum, as well as on metallicity (Gaskell & Fer- 



land| [l984 ) . The intrinsic line ratio can get as low 
as 2.6 for high metallicities, or as high as 3.4 for 
low metallicities. 

We take the metallicity into account in an ex- 
plicit manner. Starting with a fixed intrinsic 
Balmer line ratio of 3.0, we compute an initial 



metallicity estimate following Storchi-Bergmann 
et al. (1998). Once a first guess for the metal- 



licity is available, we update the intrinsic Balmer 
line ratio using an interpolation function to the 



values given in | Gaskell fc Ferland (1984). Pixel 



values now deviate individually from the initial 
guess of 3.0. With the improved extinction map 
we calculate a new abundance map. Convergence 
is achieved rapidly (less than 1% relative change 
for most pixels) after the first iteration, thus we 
stop after a second iteration. Different starting 
point values of 2.7 and 3.2 for the intrinsic Balmer 
ratio converge to the same solution. 

4.1.2. The reddening map 

The resulting reddening map is shown in the 
right panel of Fig. [4j It is based on the Seyfert 
1 mod el X-ray spectrum of Gas kell fc Ferland] 
(1984). Adopting their v~ x power law continuum 



(dashed line in their Fig. 1), E{B — V) is system- 
atically reduced by 0.16 mag, but the structures 
seen remain unchanged. We work with the higher 
extinction model for the rest of this paper, com- 
menting on the effects of a different X-ray spec- 
trum where applicable. 

The dust in J2240 is distributed unevenly. The 
area between the nucleus and the Cloud is highly 
reddened with E(B — V) = 0.5 — 0.7 mag. It is 
independent of radial velocity, and likely caused 
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Fig. 4. — Observed Ha/H/3 ratio (left), and the reddening inferred (right). 



by a large cloud of foreground dust in J2240. On 
the opposite side of the nucleus (above the dashed 
line in Fig. [4|, we find that reddening grows as 
a function of radial velocity (blue-shifted compo- 
nents are less reddened than redshifted ones). In 
this area the dust appears to be evenly embedded 
within the NLR (where it can survive, see |Laor~fc| 
Draine||1993l ), causing E(B — V) to increase from 
0.05 to ~ 0.7 mag as we look deeper into J2240. 
Other Seyfert-2 galaxies show similar (global) ex- 
tinction values of E(B — V) = 0.1 — 0.6 for their 
NLRs ( |Ho et al.|p97| |Rhee fc Larkin|[2005l ) . 



4.1.3. Error analysis 

The intrinsic Ha/H/3 ratio is affected by colli- 
sional excitation and metallicity. The latter was 
taken into account in Sect. 14.1.11 As an addi- 
tional test, we artificially increase the metallicity 
by AZ = +0.2 Z , which reduces E(B - V) by 
0.031 ± 0.005 mag only (see also Sect. [oj. The 
largest uncertainty comes from the unknown X-ray 
properties of the AGN. Using the values obtained 



by Gaskell & Ferland (1984) for a Seyfert-1 X-ray 



spectrum, we get a mean intrinsic Balmer ratio of 
2.75 with small spatial variations of 0.02 — 0.05 
mag. This is below the commonly adopted value 
of 3.08 for Seyfert-2 galaxies, and may be real, or 
the result of our ignorance about the true X-ray 
properties. The v~ x power law continuum from 



|Gaskell fc Ferland| ( p84[ ) yields 3.11±0.04 for the 
intrinsic line ratio, reducing E{B— V) by 0.12 mag, 
which is small compared to the reddening values 



up to 0.7 mag. We thus conclude that the observed 
variations in Ha/H/3 are indeed due to dust, and 
not caused by spatial variations in abundance (see 
also Sect. 4.3) nor the ionizing spectrum. 



Another source of uncertainty is the accuracy 
with which the Balmer lines have been registered. 
A Gaussian fit to the continuum profiles near Ha 
and H/3 yields an uncertainty of 0.15 pixels in the 
spatial direction. Being conservative, we shift Ha 
with respect to H/3 by 0.3 pixel in both spatial 
directions. This introduces a maximum gradual 
change in E(B — V) of 0.32 mag between +7 kpc 
and —3 kpc, about half the observed amplitude 
of 0.7 mag. Offsetting Ha by ±1 pixel along the 
spectral direction results in a 0.18 mag change be- 
tween —200 and +100 km s _1 , and no change else- 
where. This highlights the importance of an ac- 
curate registration of the Balmer lines. Assuming 
a worst-case scenario, we found that the observed 
extinction features (in particular the large dust 
cloud) prevail. We can thus realistically estimate 
a maximum registration error in E(B — V) of 0.15 
mag at any given position. 

Alternatively, E(B — V) can be obtained us- 
ing Paa/Pa/3 or other Balmer ratios. However, we 
cannot use the Paschen series as redshifted Paa 
is beyond XSHOOTER's NIR arm, and the core 
of Pa/3 is masked by airglow. Using the H7/H/3 
ratio has the advantage of not being affected by 
collisional excitation (and thus the AGN spec- 
trum). Unfortunately, redshifted H7 falls in the 
area where XSHOOTER's optical dichroic distin- 
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guishes between the UVB and VIS arms, signifi- 
cantly reducing the flux, and making flux calibra- 
tion unreliable. The S/N of the US line is too low 
to perform a 2D analysis. Integrating over the 
whole lines, both US and H7 favor an average red- 
dening of E(B - V) = 0.52 - 0.61 mag, 0.2 - 0.3 
mag higher than predicted by Ha. This can be ex- 
plained if stellar absorption, for which we do not 
correct the continuum model, reduces the flux of 
the fainter Balmer lines by 15 — 20%. 

Another consistency check is the extinction cor- 



Table 3: Electron temperatures 



rected [Om] ratio presented in Sect. 3.4 The ex- 
tinction map shows a large area below the nucleus 
of E(B — V) = 0.5 — 0.7 mag, embedded in lower 
values of 0.1 — 0.2 mag. If this differential extinc- 
tion was artificial or the [Om] ratio was left un- 
corrected, then jA5008/jA4960 increases by 0.07 for 
the highly reddened area (with respect to its sur- 
roundings). Such an increase is not found in the 
data, thus confirming this highly reddened area. 
This test is not very sensitive though, as the flux 
ratio of the relatively close [O ill] doublet is weakly 
susceptible to differential reddening effects. 

4.1.4- Testing different dust models 

The composition of dust in the harsh envi- 
ronment of an AGN is still a matter of debate, 
in particular in view of the unified dusty torus 
model (e.g. |Rhee fc Larkin 2005 and references 
therein). While dust cannot survive inside the 
broad line region (Laor & Draine 1993), larger 



grains prevail in the circum-nuclear region fur- 
ther out, leading to flat extinction curves (gray 
dust, see e.g. Maiolino et al.pOOl ). At larger dis- 
tances from the central AGN engine, smaller dust 
grains can survive or are replenished by the par- 
tial destruction of larger grains. However, very 
small grains are disintegrated quickly even at kpc 
scales (TVoit 1992 ). A large number of distant AGN 



and narrow-line quasars show extinction compati- 
ble with a dust composition like the one observed 



in the Small Magellanic Cloud (SMC; |Hopkins 
et al.||2004| |Willott|[2QQ5]) . Even steeper extinc- 



tion curves in AGN have been found by [Crenshaw 
eTaT] ( [200T| [2002] ) for NGC 3227 and Ark 564 at 
wavelengths below 4000 A. 

Applying the SMC dust model instead of R = 
3.1 causes no significant change to E(B — V) (0.04 
mag) as compared to the global variations present 
in J2240. This is expected since these models dif- 
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fer only slightly at the Ka and H/3 wavelengths 
( jCardelli et al.|1989| |Fitzpatrick|1999| ) . Both fore- 
ground dust (close to J2240) as well as dust em- 
bedded in the NLR exist in J2240. Observing lines 
further into the UV or the IR will allow us to step 
beyond the standard R = 3.1 model and constrain 
the actual physical properties of the dust in J2240, 
facilitated by the bright emission lines. 

4.2. Electron temperature and density 

The electron temperature, T e , for the medium 
ionization zone can be obtained from [Om], uti- 
lizing 



jA4960 +jA5QQ8 = 7.9 exp(3.29 X 1Q 4 /T e ) 
JA4364 1 + 4.5 xlO- 4 n e /T e - 5 



(2) 



( Osterbrock fc Ferland|2006| ) . Similar relations ex- 
ist for [Nem] ((JA3869 + jA396s)/jA3343) and [Nil] 
(C7A6548 + jA6583)/jA5755). For electron densities 
n e < 10 5 cm -3 these temperature probes become 
independent of n e . As the redshifted (and weak) 
[Nll]A5755 line is lost in the atmospheric O2 ab- 
sorption feature at 7630Awe cannot use [Nil] to 
probe the lower ionization regions. Instead, we use 
[Oh] jA3727/jA732i,32 and eq. (21) from iPeimbert 



(1967), together with an average electron density 
of 200 cm -3 for the Disk. 

Due to the weak auroral lines we cannot re- 
construct 2D temperature maps, yet we can still 
measure global values. Using [O ill] we find 17100 
K and 14700 K for the Disk and the Cloud, re- 
spectively (Table [3|. [Nem] results in 18100 K 
for the Disk, whereas the Cloud is undetected in 
[Nem]A3343 and only a lower limit of 11000 K 
can thus be obtained. Errors account for the sig- 
nificant blending of [Nelll]A3343 with the much 
brighter [Nev]A3347. For the lower ionization 
zone we find 13000 K from [Oil]. 

The [Sn] jA67i8/jA6733 ratio depends on n e and 
T e . We use T e = 13000 K as obtained from [Oil] 
for the low ionization zone. Results are shown 
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Fig. 5. — Top row: Oxygen abundance, estimated using the two strong line diagnostics of Storchi-Berg mann| 
et al.| ( p~998| ). Bottom left: Ionization parameter. Bottom right: [Sn] line ratio (map) and electron density 
(contours) for T e = 13000 K. Contour levels are drawn at n e = 150,300,450,600 cm -3 . 



in the lower right panel of Fig. [5] The Disk has 
significantly higher density than the rest, starting 
from 150 — 200 cm -3 and peaking in a dynamically 
cool spot 1 kpc north-east of the nucleus where n e 
reaches 650 cm -3 . Decreasing (increasing) T e by 
5000 K rises (lowers) the peak density by 130 (80) 
cm -3 , which is significantly less than the density 
variations measured. These features are thus real, 
and not a consequence of our simple assumption 
of constant temperature. 

The Cloud is in the low density limit of the 
[Sn] probe, thus only an upper limit of n e < 50 
cm -3 can be inferred. The area north-east of the 
nucleus, where J2240 is brightest in [S il], has lower 
density than the Disk (80 - 150 cm -3 ). 



4.3. Abundance and ionization 

Abundances for Hll regions are often obtained 



with the strong-line method (Kewley & Dopita 
2002). A more accurately approach is the direct 
method (e.g. implemented in the nebular/ionic 
task in IRAF based on the five level atomic model 



by |De Robertis et al.||1987 ). Abundance determi- 
nations for NLRs are notoriously more difficult. 



Storchi-Bergmann et al. (1998) have developed an 



empirical strong-line method based on NLRs with 
embedded Hll regions. The underlying assump- 
tion is that both the Hll regions and the NLR 
share the same metallicity, thus a calibration for 
NLRs without discernable Hll regions can be in- 
ferred. Storchi-Bergmann et al. ([1998) offer two 
methods, one based on [Nil] /Ra and [Oin]/H/3, 
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and one on log([0 ll]/[0 ill]) and log([N ll]/Ha). 
The latter depends on extinction correction as 
[O n] and [O ill] are not neighboring lines. In this 
work we also included their correction term for 
electron density. 

Both methods yield similar results of Z = 
(0.45 — 0.55) Z@, the second one with slightly larger 
values of Alog(0/H) = 0.03 (Fig.[5}. Abundances 
are enhanced by about 0.1 towards the nucleus 
and in the Cloud, and appear to increase system- 
atically with growing radial velocity. Such sub- 
solar metallicities are uncommon in NLRs, as has 



been shown by Groves et al. (2006). Their study 



of 23000 Seyfert-2 galaxies in SDSS yielded only 
~40 galaxies of evidently low metallicity. Using 
their prescription we obtain average metallicities 
of 0.9Z© (from [Nil] /[On] and [Oiii]/[Oii]), and 
- 2.0Z© (from [Nn]/Ha and [Om]/H/3). 

A large uncertainty in the abundance determi- 
nation stems from the assumptions made about 
the ionizing spectrum. Ludwig et al. (2012) re- 



port that changes to the slope and shape of the 
spectrum can change metallicities from sub-solar 
level to several times the solar value. As we do not 
know the intrinsic X-ray spectrum of the AGN in 
J2240, the true absolute value of the abundance 
remains unknown. The variations measured in the 
abundance map may be real, but can be altered 
by local variations of the X-ray spectrum filtered 
by the inhomogeneous ISM. We therefore do not 
consider abundances further in this work, leaving 
this task to future photo-ionization modeling once 
X-ray data become available. Also, note that in- 
creasing the metallicity from 0.5 Zq to several Zq 
changes the intrinsic Ha/H/3 Balmer ratio from 
2.75 to about 2.60 - 2.65, increasing E(B - V) 
only marginally. 

Once the metallicity is known, we can derive the 
ionization pa rameter U (Fig. |5| lower left) from 
[O ill] / [O ii] ( |Kewley fc Dopita||2002[) Apply ing 
the best-fit relation from |Penston et ah] ( |1990[ ) we 
find similar structures and values (not sh own). Us- 



ing [O in] /[On] requires n e < 10 3 cm 3 ( Komossa] 



fc Schulz||1997| ), since mixed-in higher densities 
will lead to an over-estimation of U. This con- 
dition is well met (Fig. [5j lower right). 

We find two strong peaks in the ionization map 
located 0.0 — 2.0 kpc north-east of the nucleus 
(above the dashed line in Fig.|5|, and spread over 
a radial velocity range of 210 km s _1 . The spa- 



tial offset is coincident with that measured for the 
density peak, and can either indicate the location 
of a second, deeply buried AGN, or simply be the 
result of a shock or interaction with a jet. Two 
weaker peaks are found 4.5 and 8 kpc south-west 
of the nucleus (below the dashed line), significant 
on the 3 and 2a level, respectively. The first is 
located between the nucleus and the Cloud in the 
area with highest dust extinction, whereas the sec- 
ond is centered in the Cloud and redshifted by 110 
km s _1 with respect to the highest ionization peak. 

The error map for U takes into account mea- 
surement errors, but not the uncertainty in metal- 
licity. To test the latter, we lower the metallic- 
ity by AZ = 0.1 Z , which corresponds to going 
from high to low metallicity areas. Log U then de- 
creases by maximally 0.12, much less than the ab- 
solute variations observed (A(logU) = 0.85), and 
less than the direct measurement errors. Note that 
while the mentioned variations in U are real, its 
absolute value should not be taken at face value, 
as the absolute metallicity is unknown. 

4.4. Continuum 

J2240 exhibits a relatively flat continuum with 
the brightest parts at 4500 A— 7000 A rest-frame 
wavelengths (Fig. 12). With an extent of about 
8 kpc in the 2D spectra, continuum radiation 
emerges from a more compact region than the line 
emission. A significant fraction of the continuum 
is stellar, as we clearly see the Call H absorption 
line at 3934A. Call K at 3969A is superimposed 
by [Nem] and He emission. The 4000A break is 
hardly visible in the data. 

No continuum can be seen in the 2D spectra 
at the position of the Cloud. However, we can 
still test for the presence of continuum flux by in- 
tegrating over all wavelengths excluding emission 
lines. The S/N obtained is sufficient to recon- 
struct broadband continuum SEDs for the Disk 
and the Cloud (Fig. [6]). While their colors are sim- 
ilar above 5500A, the Cloud is bluer than the Disk 
at shorter wavelengths. Both a nebular continuum 
and a different stellar population can cause this. 

We find a symmetric continuum profile within 
5 kpc of either side of the nucleus (Fig. [7]). In the 
north-east the profile is nearly exponential with a 
Sersic index of 0.85 ± 0.04 and a disk scale length 
of 3.0 ±0.1 kpc. In the south-west we see a long 
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Fig. 6. — By integrating the faint continuum over 
large wavelength ranges, excluding emission lines, 
we probe the photometric properties of the con- 
tinuum. The Cloud is bluer than the Disk. 

tail through the Cloud (the bump at —10 kpc in 
Fig. [7]), traceable over 18 kpc. 

4.5. Velocities and an ENLR 

The velocity FWHMs (parameter-free, directly 
measured in the line images) as a function of slit 
position for selected lines are displayed in Fig. [7| 
together with their intensity profiles. While the 
latter for low and medium ionization lines are well 
distinguished (bright and dark gray areas) within 
8 kpc of the nucleus, the velocity profiles are simi- 
lar. The notable exception is [0 11] whose FWHM 
is consistently larger by 222 ± 15 km s _1 , but oth- 
erwise has the same shape. This is a result of the 
absence of compact bright cores in [0 11] which are 
found in e.g. Ha, [Om] and [Sn] (see Fig. 
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the Appendix for the log-scaled 2D line images). 

Defining the maximum velocity as the point 
where the line profile drops to 20% of the max- 
imum emission, we find v max = (0.98 ± 0.05) x 
FWHM averaged over Ha, [Om] and [S 11] . We 
find little variation (±0.14) along the slit. Highest 
velocities are blue-shifted 1200 km s _1 . 

The [O n] velocity profile is remarkably flat be- 
tween — 14 kpc and +8 kpc, with line widths of 
534±69 km s -1 . Ha, [O ill] and [S 11] vary between 
200 and 550 km s _1 , and drop to less than 100 km 
s _1 for the largest nuclear distances. J2240 is em- 
bedded in an ENLR that co-rotates with the disk 
(albeit slower, the north-east part is redshifted 
by 195 km s _1 with respect to south-west). The 



ENLR's low luminosity, low line width, high exci- 
tation state (log([ Oin]/H/3) = 0.86-0.98) and ro - 
tation are typical flUnger et al.|1987||Du~rret|1989 ). 

Deep GMOS imaging reveals a halo stretching 
60 kpc from the nucleus (Fig. [8j top panel). To 
distinguish between a stellar tidal stream and ion- 
ized gas we use a 3h GMOS long-slit spectrum 
through a l'/5 slit, aligned along the major halo 
axis (at a position angle of 71 degrees). Unresolved 
[Oin]A5008 is detected out to 42 kpc (Fig.[8| bot- 
tom panel), thus an additional underlying stellar 
population can currently not be ruled out. The 
ENLR has a minimum diameter of 70 kpc. 

4.6. Starburst or an AGN? 



So far we have anticipated that the NLR is pow- 
ered by an AGN. Observational evidence is pre- 
sented in the following. 

Hypothetically assuming that all line emission 
is powered by stars, we can estimate the neces- 
sary star formation rate using the extrapolated 
total line fluxes from T able [4| and the scaling re- 
(119981). The result equals 



lations from Kennicutt 
290±60 M© yr _i needed to explain the Ha flux, 20 
times higher than that in an average star-forming 
GP galaxy flCardamone et al.||2009[ ). 

An AGN can be distinguished from star forma- 
tion due to its harder radiation field changing the 
line ratios. Baldwin et al. ( 1981 ) present various 



tools for this purpose. Their BPT diagrams have 
been further developed by|Veilleux fc Osterbrock 



(1987), Kewley et al. (2001), Kauffmann et al. 



(2003|) and Kewley et al. (2006). For J2240 the 



S /N is good enough to run this analysis on the ba- 
sis of individual pixels. While the precise locations 
of the dividing lines between SF, AGN, LINERs 
and composite objects in the BPT plots are still 
under some debate (|Cid Fernandes et al. 2010), 



the classification of J2240 as AGN is beyond doubt 
(Fig. [TO). [Ol]A6302,6366 emission across J2240 
shows that shock/ionization fronts are present on 
a global level, emphasizing a powerful AGN. 

When fitting the continuum across emission 
lines, we have ignored stellar absorption as the 
emission lines are much brighter than the contin- 
uum. If stellar absorption was present, Ha and 
H/3 fluxes are underestimated because of our neg- 
ligence, driving the data points in the BPT di- 
agrams away from star formation. Assuming a 
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Fig. 7. — Total flux as a function of slit position, normalized to unity, for a selection of emission lines and 
the continuum (shaded areas and thick line, respectively). The Cloud forms the bump at -10 kpc. Low 
ionization lines ([Ol], [Oil], [Nil], [Sn]) peak at 3.7 kpc from the nucleus, whereas medium ionization lines 
([Olll], [Nelll]) reach a maximum at 1 kpc. Ha runs between the two. The widths of the various lines are 
over-plotted (thin lines). While Ha, [Om] and [S ll]A6718, 33 have similar FWHMs, [Oil] is significantly 
broader. Error bars (not shown) for the velocity profiles are 10 — 15 km s _1 . 



worst-case scenario with (uncorrected) 100% stel- 
lar absorption, BPT values are still far outside the 
star-forming area (see arrows in Fig. [l0|). 

While star formation appears to be negligible 
in J2240, Bennert etH] fl2006bp warn that Hn 
regions can dominate the [O ill] flux at larger nu- 
clear distances. Levenson] ( |2007 ) and Kauffmann 



et al. (2003) also show that intense star formation 



is common in active galaxies. Indeed, a significant 
star burst can be buried by the [O ill] emission in 
J2240. However, while we do observe a decrease 
in log([Oin]/H/3) with increasing radius (middle 
left panel in Fig. [l0|), the AGN characteristics re- 
main well preserved. In addition, if the NLR is 
due to a massive outflow of hot gas, the latter may 
have disrupted star formation (for an example see 
Cano-Diaz et al.||2012[ ). 



5. A sample of GBs 

5.1. Selection and verification 

J2240 has similar colors as typical GP galax- 
ies. Applying the GP filter of |Cardamone et "aL 
(12009b to SDSS-DR8 we recover J2240 once we 



drop their maximum Petrosian radius of 2'! 0. In- 
stead, we require a minimum radius of 2'! and 
modified color cuts to select objects with particu- 
larly strong [Om] fluxes (g — r> 1.0). 95% of the 
objects found are spurious or have corrupted pho- 
tometry (the latter in particular near RA = 355±5 
and DEC = 60 ±10; see examples in Fig. 13). If in 



doubt we downloaded the original SDSS FITS im- 
ages and created the color poststamps ourselves. 
In 100% of these cases the galaxies turn out to 
have normal colors or be artifacts. 

Our SQL query for objects with 0.12 < z < 0.35 
(where [O ill] falls into the r-band) is shown in 
Appendix [Bj Shifting each bandpass by one filter 
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Fig. 8. — Deep GMOS r-band image (top) of 
J2240 revealing the extended halo. The dashed 
line indicates the position angle of the spectrum 
(bottom, showing the [Om] line), the two ar- 
rows out to which radius we detect [Om]. Black 
contour lines are logarithmically spaced. Con- 
tours (white) of the CFHT image (Fig. [I]) and the 
XSHOOTER slit (dotted line) are overplotted for 
reference. 



towards redder wavelengths (e.g. replacing g — r 
with r— i, and leaving out the constraints requiring 
filters redder than z-band), we also select possible 
candidates at 0.39 < z < 0.69. Contamination 
by lower redshift star-forming galaxies with strong 
Ha emission is expected in this higher redshift bin. 

Final manual selection results in a sample of 
29 candidates (Table [5) and Fig. [l3j. Redshifts 
and [Olll]/H/3 ratios have been measured with 
GMOS. All of our low redshift candidates observed 
show high [O m]/H/3 ratios, confirming their AGN 
character. Broad line components are absent but 
not ruled out due to shallow data. Seven galax- 
ies have a Petrosian radius between 2'! 4 and 2'/ 6 
and are good examples of GBs with particularly 
large NLRs. J2240 (listed as #016) is larger than 
any of the other galaxies (Petrosian radius of 3'/ 5). 
#004 is noteworthy, as it shows a large extension 
which indicates ejected matter or tidal interaction. 
Several GBs exhibit close companions of unknown 
redshift, indicating that encounters may play a 



role in the formation of GBs. 

In the higher redshift sample we find 12 ob- 
jects, of which 7 have spectra taken. Five galaxies 
have genuine AGN spectra, whereas 2 are star- 
forming galaxies at lower redshifts for which Ha 
was mistaken as [O ill] by our SQL filter. All ob- 
jects cover ed by the VLA FIRST catalog flWhite 
et al.|1997 ) are either non-detected or radio-quiet 

Apart from #015 all spectroscopically con- 
firmed GBs have [O ill] lines extending over 15 — 20 
kpc or more (Fig. [l4|). Given integration times of 
only 300s, the real extent of the NLRs is likely 
to be larger. In addition, the limited availability 
of guide stars is restricting the slit position angle, 
which consequently cannot be well aligned with 
the target's major axis in a majority of cases. 

5.2. Mid-IR properties 

Type-2 AGN are dusty objects with significant 
optical extinction. The WISE flWright et al.|2010[ ) 
mid-infrared fluxes can be used as a proxy for 
AGN activity, as mid-IR emission correlates with 
X-ray brightness over a wide range of luminosi- 
ties (lAsmus et al.||2011| iMason et al.||2012|). In 



particular the W4 filter at 24/im is not affected 
by dust absorption. As a comparison sample we 
choose the 887 type-2 quasars from | Reyes et al.| 
(2008), 104 of which have similar redshifts as our 
GBs (0.25 < z < 0.35). The de-redshifted mid- 
IR spectra of our GBs are very red, following a 
power-law with index (a\) = 1.99 ± 0.35, show- 
ing that the emission is likely of nuclear origin as 
compared to star-formation. The 104 type-2 com- 
parison quasars also have red spectra, albeit with 
a lower slope of (a\) — 1.59 ± 0.43. The null hy- 
pothesis that both samples have the same parent 
SED distribution is rejected on the 5% level based 
on the Kolmogorov-Smirnov test. The 24/im lu- 
minosities of the GBs and the comparison sample 
are similar though, thus GBs might simply be par- 
ticularly dusty objects. 

In Fig. |9] we plo t ^Qni] ver sus Z/24/xm for the 
obscured quasars of Reyes et al. (2008) and Greene 
( [2QTTT ). As our spectroscopic survey has not 



et al. 



been flux calibrated, we have no direct measure- 



ments of L 



[Oni] 



for the GBs. However we do 



know that the [Om] equivalent widths are com- 
parable to that of J2240, and that the spectra are 
generally similar. Assuming that [Olll]A5008 con- 



14 



45.5 



45.0 



42.5 



41 .5 



Reyes (2008) A a 

z < 0.25 . A / * * * 

• z = 0.25-0.35 Jfr A % A * ^ A 

AZ>0 ' 35 ^^fj* 

• • •••• 



^ This paper 
T Greene (201 1 ) 



45 44 45 46 

Log L [erg s~ 1 ] 

Fig. 9.— [Oiii]A5008 vs. WISE 24/im luminosi- 
ties. GBs are shown by asterisks, and the red 
triangles represent the luminous obscured quasars 
studied by Greene et al. (2011). The other data 



points show the type-2 AGN sample from |Reyes| 



et al. (2008) for three different redshift bins, with 



black dots covering the same redshift range as 
GBs. At a given L24 Mm , GBs have 5 — 50 times 
higher £[Olll] than other type-2 AGN. 

tributes the same fraction to the total r-band flux 
as for J2240 (37%), we overplot the GBs in Fig. [9] 
The 24/im luminosities for GBs and quasars in the 
same redshift range (black dots) are indistinguish- 
able. However [Om] luminosities of the GBs are 
5 — 50 times higher than expected from their mid- 
IR emission. Since the latter mainly originates 
from the compact dusty torus, this indicates that 
the current AGN activity is too low to explain the 
NLR luminosity. The NLR may therefore reflect 
an earlier, more active state, that subsided signifi- 
cantly in much less than a light crossing time. We 
return to this light echo hypothesis below. 

6. Discussion 

6.1. Main observational facts 

6.1.1. AGN character and morphology 

The galaxy-wide and ultra-luminous NLR in 
J2240 is powered by an AGN, heating the medium 
and low ionization zones to 15000 — 18000 K and 
13000 K, respectively. The Disk has higher tem- 
perature than the Cloud. The flux of the auroral 



lines is sufficiently high to facilitate more reliable 
temperature maps with deeper data. Studying the 
profile of the continuum emission within ±5 kpc 
of the nucleus we find a disk-like Sersic index of 
0.85 and a disk scale length of 3 kpc, hence ruling 
out an elliptical host galaxy. 

We find large amounts of dust as expected for 
type-2 AGN. Dust is distributed within the NLR, 
and also in a large foreground patch (Fig. [4]). Red- 
dening variations are high, and small systematics 
due to the unknown X-ray spectrum may still be 
present. BPT line diagnostics (Fig. 10) do not 
reveal significant star formation. 

J2240 is dynamically complex, showing 7 dis- 
tinct [O ill] peaks scattered over 30 kpc and 700 
km s _1 (Fig. 11). A rotational component is well 
detected in the Disk. A typical ENLR stretches 
over 26 kpc north-east and 42 kpc south-west from 
the nucleus (Fig. [8]). It can be seen in all strong 
lines such as [O n] , [O ill] , [Ne ill] , [N n] , [S n] and 
Ha, showing that J2240 is embedded in a rotating 
and quiescent, yet highly excited, bubble of gas. 

6.1.2. Shock fronts and jets 

Electron density (lower right in Fig. [5| shows 
significant fluctuations. While the Cloud is mostly 
in the low-density limit of the [S il] probe, n e rises 
sharply to 650 cm -3 in the Center. This cen- 
tral peak has a smaller FWHM (100 km s _1 ) 
than that of emission lines (200 — 500 km s _1 ), 
and it is offset by 1 kpc from the nucleus. Sev- 
eral blue- and redshifted secondary peaks exist 
within ±3 kpc of the nucleus, the redshifted ones 
forming a chain where velocity increases with nu- 
clear distance. Peaks also exist at larger separa- 
tion in the Cloud. Shock fronts, interactions with 
a jet, or mergers may cause the peaks. Indeed 
[Ol]A6302 is seen throughout the NLR, verifying 
the global presence of shocks or ionization fronts. 
Since J2240 is radio-weak, strong constraints for a 
jet are currently unavailable. The barely resolved 
radio contours in the VLA FIRST data are tenta- 
tively elongated along the major axis, thus a jet 
may explain at least some of the NLR shape. 

The ionization parameter is strongly enhanced 
in the Center and dynamically broader than the 
density. The highest ionization values are offset 
with respect to the nucleus by a similar amount 
as that observed for n e . A blue-shifted secondary 
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peak is found near the nucleus with a relative ve- 
locity of —120 km s _1 , together with two more 
ionization peaks towards the Cloud, one of which 
in the area of highest extinction. Relative veloci- 
ties are +20 and +180 km s _1 . These secondary 
peaks may be the result of shocks, or point towards 
a second AGN which may either be deeply buried 
or merely missed by our long slit observation. 

6.1.3. Globally disturbed gas kinematics 

The [O n] velocity profile is flat with FWHM = 
534 ± 69 km s _1 , and its peak coincides with the 
nucleus. Other lines, such as Ha, [O ill] and com- 
bined [S ll]A6718, 33, have almost identical profiles 
and larger variations in width than [O n] (200 — 550 
km s -1 ; Fig. [7]). The largest FWHMs of these 
lines are found 1 — 3 kpc north-east of the nucleus, 
where the medium ionization lines are bright, and 
n e and U the highest. The Cloud shows a sec- 
ondary maximum in FWHM, and a smaller one 
exists between the Cloud and the nucleus. All 
velocity profiles drop to FWHM < 50 km s _1 be- 
yond a 10 — 14 kpc radius, typical for an ENLR. 
Within that radius, maximum velocities (taken at 
20% of the line profile) are nearly the same as the 
measured FWHM. Velocities exceed 1000 km s _1 
near the nucleus, but these components are rather 
weak. Figure [7] shows that neither the gas emis- 
sion nor its velocity follow the stellar/continuum 
emission. Therefore the gas in J2240 is kinemati- 
cally disturbed out to 10 — 14 kpc radius. 

6.1.4. GB pilot survey 

Our ongoing spectroscopic pilot survey of GB 
candidates has confirmed at least 18 galaxies with 
NLRs 15 — 20 kpc in size, and those covered by the 
VLA FIRST survey are all radio-quiet. 30% of the 
GBs have close neighbors with yet unknown red- 
shifts. In several cases [Om] emission is dynam- 
ically perturbed. High-resolution spectroscopy in 
better seeing conditions than currently available 
may reveal systems with dynamic complexities as 
high as in J2240. All confirmed GBs are type- 
2 AGN, as broad-line components have not been 
identified. Upper limits to the fluxes of possible 
broad-line regions were not obtained as the spec- 
tra were taken in non-photometric conditions and 
no flux standards were observed. [O ill] luminosi- 
ties in GBs are on average one order of magnitude 



higher than expected from mid-IR emission, a sign 
of recently subsided AGN activity. 

6.2. Comparison with other NLRs 

6.2.1. Space density and [0111] luminosities 

In the 14500 deg 2 of SDSS-DR8 we find only 17 
GBs of 0.12 < z < 0.36. This redshift range corre- 
sponds to a co- moving volume of 11.1 Gpc 3 . The 
space density of GBs is thus 4.4 Gpc -3 , consistent 
with the upper end of the [O ill] luminosity func- 
tion for type-2 quasars published by Reyes ^et al.| 
fl2008p. We have neglected the 0.39 < z < 0.69 



range in this consideration due to the contamina- 
tion with star- forming galaxies. 

With 26 x 44 kpc the NLR in J2240 is much 
larger than typical NLRs (0.1 — 5 kpc, |Bennert| 



et al.||2006a| ). Approximately 99% of the [Om] 
flux are contained within a radius of 25 kpc. J2240 
falls within the scatter of the [O ill] size-luminosity 
relation observed for other Seyfert galaxies and 
quasars ( |Schmitt et aLl|2003j |Greene et aLpoTT| ). 
However, its total [Om]A5008 luminosity of 5.7 x 
10 43 ergs -1 is between 2 and 30 times higher than 
the maximum such value in the AGN samples of 



Bongiorno et al. (2010), Georgantopoulos & Aky- 
las| ( |2010D and |LaMassa et al.|d2011[ ). It also ex- 
ceeds the [O ill] luminosities of the high-z ULIRGs 
of jHarrison et al.| ([2012D by factors 2 - 12. It 



is one of the brightest objects amongst the 887 



type-2 quasars of |Reyes et al.| ([2008]). More lu- 
minous objects are only known at much higher 
redshifts, such as 2QZJ002830.4-281706 at z = 2.4 



with L^q jjjj = 2.6 x 10 44 erg s 1 (Cano-Diaz et al 



2012| ). The [Om] luminosities of the other GBs in 
our sample have so far only been estimated from 
their broad-band fluxes, but we showed they must 
also be in the range of several 10 43 erg s _1 . 

6.2.2. Radio-loud quasars and ULIRGs 

EELRs are often found around radio- loud QSOs 
(Mrk 1014 is an exception, being radio-quiet). Fu 
fc Stockton| (12009 ) identify compact clouds mov- 
ing with high velocities (~ 500 km s _1 ), yet having 
low line widths (30 — 100 km s _1 ). They do not 
show morphological links with their host galax- 
ies. EELRs are probably formed during mergers, 
which start the QSO engine that ionizes the gas 
and blasts it into the outer surroundings. 
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J2240 is different. Firstly, its host galaxy is 
radio- weak or -quiet. Secondly, the NLR is at least 
one order of magnitude brighter than the EELRs 



described in |Fu et al-| ( [2012| >. Thirdly, J2240 shows 
high line widths (200 — 500 km s _1 ) and low ra- 
dial velocities (~ 200 km s _1 ), opposite of what we 
see in EELRs. This may be an effect of incomplete 
sampling by the long slit, as we might have missed 
the high velocity components. Or the gas clouds 
move predominantly perpendicular to the line of 
sight such that we do not see the full velocity vec- 
tor. Lastly, the NLR embeds the host galaxy, thus 
establishing a morphological link. The same dif- 
ferences are also found when comparing GBs to 
the galaxy-scale emission line regions observed in 
powerful radio galaxies and some ULIRGs in the 
young Universe (z = 2 — 3; Nesvadba et al.|[2008 



Harrison et al. 2012). This holds for the other GBs 



as well given our current data. 

6.2.3. Current light echo samples 

J2240 does not appear to be a typical light 
echo, such as observed in Hanny's Voorwerp (Lin- 
tott et aL]|2009 ), or those described by Keel et al. 
(2012a). There are similarities though. Firstly, 



our BPT diagrams (Fig. 10) are fully compatible 
with the ones of Keel et al. (2012a). Secondly, in- 



teraction with neighboring galaxies are frequently 
seen in the sample of Keel et al. (2012a), and com- 



panions are also observed for the GBs. 

Nevertheless, typical [O ill] luminosities for 
these light echos are 2 orders of magnitude lower 
than for GBs. This is also evident in the SDSS im- 
ages, which show the stellar body of these galaxies 
well, whereas for GBs the stellar emission is over- 



whelmed by the NLR. The survey of Keel et al. 
(2012a) gets insensitive at redshifts z > 0.1, as 
the investigated features become too faint. Our 
GB sample on the other hand, extracted from the 
same data base, does not reveal any GBs with 
z < 0.19 (we are sensitive down to z = 0.12). 
Shifting our SQL filter to even lower redshifts did 
not reveal further candidates. This is consistent 
for two reasons. First, a galaxy like J2240 at 
redshifts less than 0.1 would be conspicuous and 
likely picked up by earlier surveys. Second, the 
comoving volume within z < 0.1 is just about 
0.3 Gpc 3 . If the space density of GBs calculated 
above does not evolve between z ~ 0.3 and today, 
then we expect only 0.44 GBs at z < 0.1 within 



the 14500 deg 2 cove r ed by S DSS-DR8. For com- 
parison, Keel et al. (2012a) find ~ 100 galaxies 



with possible light echos at z < 0.1. 

6.3. Binary AGN or SMBH merger? 

Can the extraordinary properties of J2240 be 
explained by mergers? We observe globally dis- 
turbed gas kinematics, and the highly ionized com- 
pact [O ill] sources can indicate multiplicity (|Com-| 
erford et al. 2012). A tidally distorted neighbor- 
ing galaxy (Fig. |1| makes a multiple merger sce- 
nario plausible, in which another, possibly gas- 
rich, galaxy is currently consumed by J2240. As 
Liu et al. (|2012|) demonstrate for wide separa- 



tion binary AGN, the SMBH accretion rate is in- 
creased in such double systems as the merger pro- 
cess funnels more material towards the centers. 
LogQOlll]) increases by 0.7 ± 0.1 when decreas- 
ing the separation from 100 to 5 kpc in these sys- 
tems. Even closer pairs likely have correspond- 
ingly higher [O ill] luminosities. 

Are we witnessing some violent process during 
the final stages of an AGN merger? SMBHs are 
common in the centers of massive galaxies, thus 
galaxy mergers must also result in the coalescence 
of SMBHs. Accordingly, SMBH or AGN pairs 
should be common. However, with decreasing sep- 
aration they are increasingly hard to identify. For 
example, numerous binary AGN with separations 
of tens of kpc are known, yet with 3.6% their 
fraction among optically selected AGN is already 
small (|Liu et al.l|2011[) |Shen et al.| ( |2011[ ) and 
Comerford et al. fl2012|) find kpc binary AGN in 



galaxies with double-peaked [O ill] emission, how- 
ever that feature is more commonly caused by gas 
kinematics and cannot be used as a reliable indi- 
cator for AGN binarity ( |Shen et al.|2011| |Fu et al. 
|2012|). X-ray, infrared or radio observations are 



needed to confirm such systems. 

Only few binary AGN with even smaller separa- 



tions are known ( Rodriguez et al.||2006 Fabbiano 
et al.pOlT ), including a system with sub-pc scale 
and an orbital period of ~ 100 years (Boroson & 



Lauer 



2009) 



Eracleous et al. 



(2012) report sev- 



eral sub-pc candidates, however emphasizing that 
long-term monitoring is required to link observed 
line variability with orbital motions. From the X- 
ray perspective statistics are equally weak. For 
example, |Teng et al.| ( [2012] ) find that only - 8% 



of massive mergers actually harbor binary AGN. 
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While the observational data base is poor, simu- 
lations of the pre-coalescence state of SMBH merg- 
ers have been carried out (Hopkins & Quataert 
20101 | Khan et ah] [2012] |Van Wassenhove et a" 



2012). During the orbital decay from kpc to pc 
scales no processes are found that explain an AGN 
flaring up for 10 4 — 10 5 years by 3 — 4 orders of 
magnitude. Neither can current accretion models 
explain a shut-down on similarly short time scales 
dSchawinski et al.|201Q ). 

It is also unlikely that we are observing some ef- 
fect or aftermath of the actual coalescence of two 



SMBHs. Tanaka et al. (2010) calculate the electro- 



magnetic footprint of SMBH mergers. They find 
an increase of bolometric luminosity of about 10% 
per year over time scales of years or decades, to- 
gether with an increase of X-ray hardness. This 
is, however, much less than the light crossing time 
of a galaxy (10 4 — 10 5 years) needed to explain the 
size and luminosity of our NLRs. 

6.4. Light echos - quasars shutting down? 

The most prominent characteristics of the 
NLRs in our GBs are their large angular extent 
and high [Om] luminosities. Certainly powerful 
AGN must be responsible for this, but they are 
not evident from the SDSS imaging. Either they 
are deeply buried, or their activity has steeply de- 
clined over time scales much less than the light 
crossing time of the NLR, in which case we are 
observing strong light echos. 

To this end we have to show that the current 
AGN activity is much lower than expected from 
the overall [Om] luminosity, and that we are not 
just observing very obscured nuclei. The best way 
of doing this is to determine the current X-ray 
luminosities of GBs, and compare them to their 
[Om] luminosity. In case of a light echo, the 
X-ray output will not match the [Om] luminos- 
ity. Using the best-fit relation between Lx and 
£[Oni] f° r type-2 quasars and Seyferts, we expect 



Lx ~ lx 10 44 erg s 1 in the 2 — 10 keV range Oia 



et al.|2012||LaMassa et al.|2009[ ). Note that, since 
ROSAT is only sensitive to soft energies, even the 
lowest plausible X-ray absorption is sufficient to 
account for the non-detection of J2240 by ROSAT. 
Therefore we cannot constrain a possible quasar 
shutdown with the X-ray data currently available. 
To overcome the lack of X-ray data we use the 



mid-IR emission as a proxy for AGN activity, as 
it is unaffected by dust absorption and emanates 
from the immediate, pc-scale AGN environment. 
Mid-IR emission is tightly correlated with the X- 
ray luminosity in AGN ( |Asmus et al.|20lT Mason 



et al.||2012| ). We compare WISE 24/im luminosi- 



ties with [Olll]A5008 luminosities for a large sam- 
ple of type-2 quasars and GBs. Any change in 
AGN activity in GBs will need about a galaxy's 
light-crossing time before it is fully reflected in the 
NLRs' properties. We find the [O ill] luminosities 
of the GBs to be 5 — 50 times higher than expected 
from the mid-IR emission, strengthening the light 
echo scenario. If J2240 and the other GBs are in- 
deed light echos, they are spectacular examples of 
powerful QSOs currently shutting down. 

Note, however, that the [Om] luminosities of 
our GBs surpass the highest [O ill] luminosities of 
the reference type-2 sample at the same redshift 
(Fig. [9]). If the light echo interpretation is correct, 
then one expects quasars at similar redshifts with 
suitable mid-IR fluxes that have not shut down 
yet. Only quasars at z > 0.4 match such high 
fluxes. We think that this is a selection effect, 
with GBs dropping out of the main SDSS spectro- 
scopic target selection algorithm. The fact that 
our simple size and color selection of GBs yields 
about 95% spurious sources shows that the combi- 
nation of broad-band photometry and angular size 
is quite unusual. The SDSS algorithm designed to 
select sources for spectroscopic follow-up is opti- 
mized for a high success rate, allowing for only 
5% unusual sources to be included (see e.g. |Reyes| 
et al. 2008, and references therein). It is there- 
fore not surprising that only one of our GBs has a 
SDSS-DR8 spectrum. 

6.5. Outlook 

Considering the intrinsic luminosity that must 
be responsible for the observed optical NLR, the 
hard X-ray continuum of AGN must be directly 
detectable, provided that the column density is 
not too great (Nn < 10 23 cm -2 ). Even if the cen- 
tral engine is deeply obscured or truly hidden by 
Compton-thick absorption, the characteristic Fe 
Ka fluorescence line should be detectable. In this 
way a binary AGN with a few kpc separation can 
be directly confirmed through X-ray imaging, or 
the shut-down time-scale of the obscured quasars 
in GBs constrained. X-ray observations will also 
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determine the slope of the ionizing spectrum, im- 
proving photo-ionization models of these NLRs. 

Additional work going beyond our simple long- 
slit spectroscopic analysis has to be done. A real- 
istic model of the NLR requires the use of photo- 
ionization codes and knowledge about the X-ray 
properties. The optical spectral analysis must also 
be extended to the full body of J2240. Such in- 
vestigation may reveal ionizing sources missed by 
the long-slit observation, and provide better ex- 
tinction, temperature and density maps. In addi- 
tion we can better constrain the dynamics, deter- 
mine the gas mass, and investigate whether GBs 
have significantly higher gas masses than other ob- 
scured quasars showing massive outflows or large 
NLRs. To this end, GMOS IFU observations in 
good seeing conditions have been conducted of 
J2240 in 2012B at the Gemini Observatory. We 
have also applied for follow-up observations of sev- 
eral more GBs. 
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Further J2240-0927 data 



Table 4: Fluxes (corrected for galactic extinction) and luminosities of J2240 for selected emission lines, and 
their relative strength with respect to H/3. Total values are corrected for slit losses by seeing and for the fact 
that the slit is smaller than the object. The correction factor used is 1.9 ±0.3, its uncertainty conservatively 
estimated. 



Line 


jslit 

[xlO -16 erg s _1 cm -2 ] 


L slit 

[xlO 40 erg s- 1 ] 


L tot 

[xlO 40 erg s- 1 ] 




Ha 


338.6 ±2.8 


1200.2 ±9.9 


2280 ± 360 


4.25 


H/3 


79.7 ± 1.2 


282.4 ±4.1 


536 ± 85 


1.00 


H 7 


31.6 ±0.9 


111.9 ±3.3 


212 ± 34 


0.40 


RS 


12.8 ±0.3 


45.6 ± 1.1 


86 ±13 


0.16 


HeiA5877 


7.6 ±0.1 


26.7 ±0.6 


50 ±8 


0.09 


HeiiA4687 


10.6 ±0.4 


37.5 ±1.3 


70 ±11 


0.13 


HeiA10833 


23.7 ±0.8 


84.0 ±2.9 


159 ± 25 


0.30 


[0 1] A6302 


44.0 ±0.6 


156.0 ±2.1 


296 ± 46 


0.55 


[0 n] A3727 


296.6 ±4.1 


1051.3 ±14.7 


1997 ±316 


3.72 


[0 m] A4364 


10.3 ±0.5 


36.4 ±1.9 


69 ±11 


0.13 


[0 in] A5008 


841.5 ±14.0 


2982.7 ±49.7 


5667 ± 899 


10.56 


[Nil] A6586 


187.7 ±1.6 


665.3 ±5.5 


1264 ± 199 


2.36 


[S n] A7618,33 


168.1 ± 1.6 


595.6 ±5.4 


1131 ± 178 


2.11 


[S in] A9071 


29.6 ±0.9 


104.9 ±3.1 


199 ± 32 


0.37 


[S in] A9533 


107.2 ±2.6 


379.8 ±9.2 


721 ± 115 


1.34 


[Nem] A3869 


56.3 ±0.9 


199.6 ±3.0 


379 ± 60 


0.71 


[Nev] A3427 


17.4 ±0.4 


61.8 ±1.6 


117 ± 18 


0.22 


[Mgn] A2799 


19.5 ±0.5 


69.1 ± 1.9 


131 ±21 


0.24 


[Arm] A7138 


8.9 ±0.1 


31.4 ±0.6 


59 ±9 


0.11 
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Fig. 10. — BPT diagrams. Top: The solid and dotted lines represent the model division between the two 

). Data points show individual pixels in the NLR, their 
mean errors are indicated. The NLR is evidently powered by an AGN. The maximum possible effect of 
stellar absorption is indicated by the arrows. Middle and bottom: Spatially and spectrally resolved BPT 
diagrams, based on the same data points as shown in the top row. 



galaxy types and its uncertainty ( Kewley et al.pOOl 
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Fig. 11. — Top 4 rows: Selected emission lines (continuum subtracted, and projected to a common reference 
frame) in J2240. A logarithmic scaling (different between lines) has been chosen for better visualization. 
The dashed line indicates the spatial position of the nucleus (position of highest continuum level), i.e. the 
spatial direction is oriented along the vertical axis. Wavelength increases from left to right. Spatial and 
velocity scales are indicated in the upper left panel. The area displayed for each line is slightly larger than 
that used for the 2D spectral analysis plots, in order to show more details. Bottom: For comparison, we 
show a linearly scaled larger part of the spectrum before continuum subtraction, centered on [Oll]A3727. 
This gives an idea about the strength and extent of the emission lines (see also Fig. 12). 
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Fig. 12. — Rest-frame ID XSHOOTER spectra of J2240, uncorrected for galactic extinction. The black line 
represents the galaxy center, integrated within ±4.5 kpc of the nucleus. The blue line has been integrated 
over 7.6 kpc centered on the Cloud. Note the great similarity between the two spectra. For visualization 
purposes data have been filtered with a 7 A wide median kernel, thus the actual resolution is 48 (12) times 
higher than shown for the UVB/VIS (NIR) channels^espectively. 



B. The GB sample 

The following SQL filter was used to retrieve the initial selection of 376 GB candidates from SDSS- 
DR^] including poststamps for a visual cross-check. Only ~ 5% of the objects are retained as genuine GB 
candidates after visual inspection. We experimented with various parameters, such as the Clean flag, but 
found that these were prone to exclude genuine objects (amongst others J2240). 



SELECT 



FROM 



WHERE 



ra, dec, objID, u,g,r,i,z, petrorad_r, 

; <a href =http : //cas . sdss . org/dr3/en/tools/chart/navi . asp?ra= ; + 

cast(ra as varchar(10))+ , &dec= , +cast( dec as varchar(lO)) + >>> + 

' <img src="http : / / skyservice . pha . jhu . edu/dr8/ImgCutout/get jpeg . aspx?ra= } + 

cast(ra as varchar(15))+ , &dec= , +cast(dec as varchar (15) )+ 

, &scale=0.40&width=120&height=120&opt="/> ' as pic 



Galaxy 




((r 


<= 20.5) 




and 


(r >= 17) 




and 


(u - r <= 5) 




and 


(r - i <= -0.2) 




and 


(r - z <= 0.6) 




and 


(g - r >= r - i 


+ 0.5) 


and 


(u - r >= 2.5 * 


(r - z)) 


and 


(g - r > 1.0) 




and 


(petrorad_r > 2) 




and 


(psfmagerr_g < 


.04) 


and 


(psfmagerr_r < 


.04) 


and 


(psfmagerr_i < 


.04)) 



t http: / / skyserver.sdss3.org/dr8 /en /tools /search /sql.asp 
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Fig. 13.— Top three rows: GBs with 0.12 < z < 0.36 (see also Table [5}. #016 is J2240. Rows 4 and 5: 
Possible GBs with 0.39 < z < 0.69. Bottom: common false positives, from left to right: bad background, 
systematic photometry error or bad seeing, bright star/diffraction spike, Herbig-Haro object, satellite, star 
near detector edge. 
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Fig. 14. — [Olll]A5008 (log-scaled) for some GBs. The continuum has not been subtracted. 
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